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Rough Gauge Fields, Smearing and Domain Wall Fermions 

Meifeng Lin a [RBC Collaboration] * 
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At a fixed lattice spacing, as determined by say m p , adding additional fermion flavors to a dynamical simu- 
lation produces rougher gauge field configurations at the lattice scale. For domain wall fermions, these rough 
configurations lead to larger residual chiral symmetry breaking and larger values for the residual masses, m rcs . 
We discuss ongoing attempts to reduce chiral symmetry breaking for Nf = 3 dynamical domain wall fermion 
simulations by different smoothing choices for the gauge fields. 
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1. INTRODUCTION 

As is widely discussed (see, for example, 
Rcf. 1 ), topological lattice dislocations give rise 
to a non-vanishing density of near-zero modes of 
the hermitian Wilson-Dirac operator (HWDO), 
which is considered to be an origin of large chi- 
ral symmetry breaking for domain wall fermions 
(DWF) and a sign of the theory being close to 
or inside the Aoki phase. In quenched simula- 
tions, we are able to suppress the creation of these 
lattice dislocations by using improved gauge ac- 
tions 2 (DBW2 for instance) and have achieved 
better chiral symmetry in numerical simulations. 
However, if the lattice spacing, as determined by, 
say, the p mass, is kept fixed as more flavors of 
fermions are added, asymptotic freedom predicts 
a rougher gauge field at the single-lattice spacing 
scale. As a consequence of the rougher fields, the 
residual mass m res 3. , becomes larger at a compa- 
rable lattice spacing than in the quenched theory. 

Recently there have been various attempts ^ 
to improve the chiral properties of DWF in dy- 
namical simulations by modifying the gauge ac- 
tion. Here we study a smearing technique on 
Nf — 3 lattices, i.e., replacing each gauge link 
by a weighted sum of the original link and prod- 
ucts of gauge matrices along 3, 5 and 7 link paths, 
to see if chiral symmetry breaking is reduced. 



* We thank RIKEN, Brookhaven National Laboratory and 
the U.S. Department of Energy for providing the facilities 
essential for the completion of this work. 



2. SMEARING ON GAUGE LINKS 

A smeared link is constructed as: 



+ c 3 £ 



+ c 5 E 




+ c 7 Z 




As we can see, the new links are not necessarily 
elements of £17(3). Although this may make the 
interpretation of the unrenormalized observables 
obscure, it does have the advantages of cheap 
computer time cost and simple numerical imple- 
mentation in full QCD simulations. Yet the c's 
should be appropriately chosen so that the con- 
ventional field normalization in weak coupling is 
preserved. 

The c's can be simply normalized as: 



ci + 6c 3 + 24c 5 + 48c 7 = 1 



(1) 



While this preserves the field normalization in 
weak coupling, at the lattice spacings of current 
simulations a better choice has been seen to be 
the so-called tadpole improvement condition, in 
which the smearing coefficients are chosen so that 



c\ + 6c 3 u n 



24c 5 Uo 



48e 7 i^ = 1 



(2) 



where u$ is the quartic root of the average pla- 
quette. We explore these and other normaliza- 
tion possibilities by measuring m res with a fixed 
M5 = 1.8 for different combinations of c\ and 
C3 without enforcing a pre-defined normalization. 
Coefficients which give smaller m res are then fur- 
ther studied. 
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Table 1 

m res with different smearings. All the measure- 
ments have valence L s — 8 and M 5 = 1.8 except 
the one with c-j ^ 0, which has M 5 = 2.2. Data 
for unsmeared lattices are from 84 configurations; 
others are all from 30 configurations. 
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3. SIMULATION RESULTS 

3.1. Residual masses 

Smearings are performed on our 16 3 x 32, 
Nf — 3 lattices, which were generated with the 
DWF and DBW2 gauge action. The parame- 
ters used were (5 = 0.72, M 5 = 1.8, L s = 8 
and r7j soa = 0.04, yielding a lattice spacing (from 
m p ) of a- 1 w 1.6 to 1.7 GeV 0. Choices 
for the smearing coefficients we have studied in- 
clude (i) one set normalized according to eq.Q: 
ci = g,c 3 = ^,c 5 = ^,c 7 = 3^,0 la MILC's 
Fat7; (ii) one set following eq.@: c\ — 0.25, c 5 = 
0.051; (iii)&(iv) two sets with arbitrary normal- 
izations for comparison: ci — 0.4, c 3 = 0.12 and 
ci = 0.8, c 3 = 0.06. In Table □ we see that 
TO res has generally decreased for different choices 
of smearing coefficients except set (i) . 

For DWF, the value of 

777-res IS affected by the 
value for M5. However, smearing can change the 
normalization of terms in the action which involve 
the gauge field, relative to the normalization of 
the mass term M 5 . We can see this by consid- 
ering eigenspectra of the HWDO (spectral flows) 
on the smeared gauge configurations. Shown in 
Figure Q are representative spectral flows from 
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Figure 1. Eigenspectra of the hermitian Wilson- 
Dirac operator for different smearing coefficients 
on the same lattice. 



smeared and unsmeared lattices 1 . The pinch 
point (k c for Wilson fermions) on the top graph 
(Fat7 smearing) has shifted from — M5 ~ —0.9 
to ~ —1.8. The larger m ros for this smearing 
choice is a result of non-optimal M5 and possi- 
bly rougher gauge configurations. On the bottom 
graph (tadpole smearing), the density of small 
eigenvalues is visibly decreased and the gap is 
widely opened near — M5 = —1.8, leading to the 
smaller residual masses with coefficient choice (ii). 
However, it is also worth noting that the major 
crossing at —A/5 ~ —1.6 still persists after smear- 
ing. The failure to eliminate these crossings indi- 
cates that the topological lattice dislocations are 
not changed. While m res has decreased for this 
L s , the persistence of the topological dislocations 
makes studying the large-L s behavior of m res im- 
portant. Wc present results on this in Sec. 13.31 



lr rhe author thanks S.D. Cohen for providing the data nec- 
essary for this figure. 
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Table 2 

m p a in valence quark chiral limit and mass 
renormalization. m va i = 0.02, 0.03, 0.04, 0.05 are 
used in the extrapolation for unsmeared lattices; 
m vai = 0.02, 0.04 are used otherwise. 

Cl Cj C 5 C7 TO p a Z m 

1.0 0.539(4) 1.0 

g 4 B ^ 0-522(28) 0.76(3) 

0.25 0.051 0.564(34) 0.79(2) 

0.4 0.12 0.563(31) 0.88(2) 

0.8 0.06 0.563(26) 0.83(3) 
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Figure 2. L s dependence of m res . 

ual masses show simple exponential decay with 
L Sl here both smeared and unsmeared lattices 
show a more complicated behavior. The same 
phenomenon has also been seen in our quenched 
simulations with Wilson gauge action pj], where 
many topological dislocations are also observed. 

4. CONCLUSIONS 

We have found that for DWF, the smearings we 
have chosen do not change the lattice scale, but 
introduce additional mass renormalization fac- 
tors. The smeared gauge fields are smoother, but 
we conclude that the smearings we have chosen 
do not change the topological dislocations present 
on the dynamical lattices, thus the large L s be- 
havior of m rcs stays the same. Smearing does 
not appear to help DWF simulations at strong 
coupling. However, for a^ 1 ~ 2 GeV, where dis- 
locations are suppressed due to weaker coupling, 
and a practical choice of L s (say, 16), dynamical 
3 flavor simulations with DWF are ready to do 
with current techniques. 
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3.2. Lattice scales and mass renormaliza- 
tion 

To solidify our conclusions about the effects 
of smearing on m rcs , we have studied the lattice 
spacing on smeared lattices using m p , since it is 
known that m res is strongly dependent on lattice 
spacing. We determine the lattice scales by tak- 
ing m p a in the valence quark chiral limit. As 
shown in Table |3 m p a is essentially unaffected 
by smearings. 

Another important effect when comparing m res 
with different smearing choices is the change in 
mass renormalization introduced by smearing. 
Noting the fact that m^. ex m 9 , where m q = 
TO va i + m res for DWF, we can extract the smeared 
lattice mass renormalization factors relative to 
those on the unsmeared lattices by 

(mt) 2 lm s „ 

Z ™ = "V q (3) 

where the quantities on unsmeared lattices are 
superscripted by u and those on smeared lattices 
by s. The renormalizcd residual masses are thus 
m£g S = Z m m ves . The Z m 's are also shown in Ta- 
ble 

3.3. L s dependence of m rcs 

Given that smearing has not removed the cross- 
ings at —M$ = —1.6 (lower panel of Fig. 
we have done further investigations on the de- 
cay of m rcs with L s . As shown in Figure [21 
m res on the smeared lattices decays essentially 
the same as on the unsmeared lattices. Contrary 
to the quenched DBW2 case [2, where the resid- 
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